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Abstract 

This  paper  examines  the  realization  of  acausal  weighting  patterns  with  two-point 
boundary-value  descriptor  systems  (TPBVDSs).  We  restrict  our  attention  to  the 
subclass  of  TPBVDSs  which  are  extendible,  i.e.,  whose  input-output  weighting 
pattern  can  be  extended  outwards  indefinitely,  and  stationary,  so  that  their 
weighting  pattern  is  shift-invariant.  Then,  given  an  infinite  acausal  shift- invariant 
weighting  pattern,  the  realization  problem  consists  in  constructing  a  minimal 
TPBVDS  over  a  fixed  interval,  whose  extended  weighting  pattern  matches  the 
given  pattern.  The  realization  method  which  is  proposed  relies  on  a  new 
transform,  the  (s,i)  transform,  which  is  used  to  determine  the  dimension  of  a 
minimal  realization,  and  to  construct  a  minimal  realization  by  factoring  two 
homogeneous  rational  matrices  in  the  variables  s  and  t . 

1.  Introduction 

There  exists  an  extensive  literature  [l]-[3]  on  the  state-space  realization  prob¬ 
lem  for  linear  time-invariant  causal  systems,  i.e.,  for  systems  which  admit  an 
input-output  description  of  the  form 

OO 

y{k)=  E  W(A:-/)u(0  (1.1) 

/=— OO 

where  the  impulse  response  (weighting  pattern)  W (.)  satisfies 

W{m)  —  0  for  m  <  0  .  (1-2) 

However,  for  many  physical  systems,  in  particular  when  the  independent  variable 
is  space  rather  than  time,  the  causality  condition  (1.2)  does  not  hold.  For  exam¬ 
ple,  if  we  consider  the  temperature  of  a  heated  rod,  there  is  no  reason  to  assume 
that  the  temperature  at  any  point  of  the  rod  depends  exclusively  on  the  applied 
heat  on  one  side  of  that  point.  Weighting  patterns  that  do  not  satisfy  (1.2)  are 
called  acausal.  In  this  paper,  we  develop  a  realization  theory  for  acausal  weighting 
patterns  in  terms  of  two-point  boundary-value  descriptor  systems  (TPBVDSs). 

The  motivation  for  considering  this  class  of  systems  is  that  the  dynamics  of 
discrete-time  descriptor  systems  are  noncausal,  in  the  sense  that  they  contain  com¬ 
ponents  which  propagate  in  both  time  directions  [4].  The  boundary  conditions  are 
another  source  of  noncausality,  since  they  are  expressed  symmetrically  in  terms  of 
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the  system  variables  at  both  ends  of  the  interval  of  definition.  Thus,  TPBVDSs 
have  a  totally  acausal  structure  [5]- [6],  which  is  ideally  suited  to  model  noncausal 
systems.  Motivated  by  the  earlier  work  of  Krener  [7],  and  Gohberg,  Kaashoek  and 
Lerer  [8]  for  acausal  systems  with  standard  dynamics,  a  complete  system  theory  of 
TPBVDSs  has  been  developed  recently  in  [9]-[ll],  including  concepts  such  as 
reachability,  observability  and  minimality.  In  this  paper,  we  restrict  our  attention 
to  stationary  and  extendible  TPBVDSs,  i.e,  TPBVDSs  whose  weighting  pattern 
pattern  is  shift-invariant,  and  where  the  interval  of  definition  of  the  TPBVDS  can 
be  extended  outwards  indefinitely,  without  changing  the  weighting  pattern. 

The  realization  problem  that  we  consider  can  be  stated  as  follows;  given  an 
infinite  weighting  pattern  W{m),  construct  a  minimal  TPBVDS  over  a  fixed  inter¬ 
val,  which  has  W{m)  for  extended  pattern.  As  for  causal  time-invariant  systems, 
where  the  .z-transform  plays  a  useful  role  in  transforming  the  realization  problem 
into  a  factorization  problem  for  proper  rational  matrix  transfer  functions,  it  is 
shown  that  the  TPBVDS  realization  problem  can  be  formulated  as  a  factorization 
problem.  However,  instead  of  using  the  .^-transform,  we  introduce  a  new 
transform,  the  {s,i)  transform,  which  handles  zero  and  and  infinite  frequencies 
symmetrically,  and  is  therefore  well  adapted  to  the  analysis  of  descriptor  systems. 
The  {s,t)  transform  is  used  here  to  characterize  the  dimension  of  TPBVDS  reali¬ 
zations  in  terms  of  the  McMillan  degree  for  rational  matrices  in  s  and  t,  and  to 
formulate  the  TPBVDS  realization  problem  as  a  factorization  problem  for  not 
one,  but  two  homogeneous  rational  matrices  in  two  variables.  Due  to  space  limita¬ 
tions,  most  results  are  stated  without  proof.  The  reader  is  referred  to  [12]  for  a 
complete  account. 

2.  Two-point  Boundary- value  Descriptor  Systems 

In  this  section,  we  review  several  properties  of  TPBVDSs,  such  as  stationar- 
ity,  minimality  and  extendibility,  that  will  be  needed  in  the  development  of  our 
TPBVDS  realization  procedure. 

2.1.  Model  Description 

A  linear  time-invariant  TPBVDS  is  described  by  the  difference  equation 

Ex{k+1)  =  Ax{k)  +  Bu{k)  ,  0<k<N-l  (2.1) 

with  boundary  condition 

y,  x(0)  -t-  Vyx(iV)  =  t;  (2.2) 

and  output  equation 

y{k)  =  Cx{k),  0<k<N.  (2.3) 

Here,  x  and  v  are  n-dimensional,  u  is  m-dimensional,  y  is  p-dimensional,  and  E , 
A,  B  and  C  are  constant  matrices.  We  also  assume  that  N  '>2n,  so  that  all 
modes  can  be  excited  and  observed.  In  [9]  it  was  shown  that  if  the  system  (2.1)- 
(2.2)  is  well-posed,  we  can  assume  without  loss  of  generality  that  (2.l)-(2.2)  is  in 
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normalized  form,  i.e.,  that  there  exists  scalars  a  and  /?  such  that 

OiE  -f-  ^A.  =  I  (2*4) 

(this  is  referred  to  as  the  standard  form  for  the  pencil  {E ,  yl })  and  in  addition 

ViE^  .  (2.5) 


Then,  the  map  from  {u,t;}  to  x  has  the  form 

x{k)  =  A^ E^-’^  V  +  G{k,l)Bu{l)  ,  (2.6) 

where  G{k,l)  is  the  Green’s  function  associated  to  the  descriptor  dynamics  (2.1) 
and  boundary  condition  (2.2).  The  map  from  inputs  u  to  outputs  y  specifies  the 
weighting  pattern  W  of  the  system.  Setting  u  =  0  in  (2.6)  yields 

y{k)~'EW{k,l)„(l),  (2.7) 

1=0 


with 


W{k,l)  =  CG{k,l)B  . 


(2.8) 


2.2.  Stationarity 

In  contrast  with  the  causal  case,  where  time-invariant  state-space  models 
have  a  time- invariant  impulse  response,  the  weighting  pattern  W{k,l)  is  not  in 
general  a  function  of  the  difi'erence  k—l.  TPBVDSs  that  have  this  property  are 
called  stationary. 

Theorem  2.1  [10]:  The  TPBVDS  (2.l)-(2.3)  is  stationary  if  and  only  if 


0,[Vt,E\R,^0,\V,^A]R, 

(2.9a) 

0.[V,,E]R,=0,\Vi_A]R,  , 

(2.9b) 

where  [X, 

T]  denotes  the  commutator  product  of  X  and  Y 

[X,  Y]=XY  -  YX 

(2.10) 

and 

- 

R,  =  [E’^-^B  AE’^-^B  •  •  ■  A’^-^B] 

(2.11a) 

0/=  [{E^-YC^  {AE^-^)^C'^  •  •  •  (A"-^)^^^]  . 

(2.11b) 

The  matrices  and  in  (2.11)  are  respectively  the  strong  reachability  and 
strong  observability  matrices  of  the  TPBVDS,  as  discussed  in  [9].  The  stationarity 
conditions  (2.9a)  and  (2.9b)  state  that  V,-  and  must  commute  with  E  and  A , 
except  for  parts  that  are  either  in  the  left  null  space  of  or  the  right  null  space 
of  Og.  Consequently,  if  R^  and  0^  have  full  rank,  i.e.,  if  the  TPBVDS  is  strongly 
reachable  and  strongly  observable,  F,-  and  Vj  must  commute  with  E  and  A  . 

It  is  shown  in  [10]  that  the  weighting  pattern  of  a  stationary  TPBADS 
defined  over  [0,V]  is  given  by 
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CViA’^-^E 

CVfE~’^A 


N-kB 
N +k  —  \ B 


l<k  <N 
1-N  <k  <0  . 


(2.12) 


2.3.  Minimality 

Since  our  goal  is  to  realize  shift-invariant  acausal  weighting  patterns  with 
stationary  TPBVDSs,  we  need  to  be  able  to  determine  whether  a  system  in  this 
class  is  minimal  or  not.  This  issue  was  studied  in  detail  in  [10]-[11],  leading  to  the 
following  characterization  of  minimality. 

Theorem  2.2:  The  stationary  TPBVDS  (2.l)-(2.3)  is  minimal  if  and  only  if 


(a) 

(b) 


[V^Rg  VjRg]  has  full  row  rank, 

fo.nl 


O.V, 


f  J 


has  full  column  rank, 


(c)  Ker(OjClm(i2j- 


(2.13a) 

(2.13b) 

(2.13c) 


It  was  also  shown  in  [10,  Corollary  5.1]  that  Theorem  2.2  implies: 

Corollary:  Let  ,Ej,Aj-,N)  with  j  =  1,2  be  two  minimal  and  sta¬ 

tionary  realizations  of  the  same  weighting  pattern,  where  {Ej,  Ay},  j  —  1,2  are  in 
standard  form  for  the  same  a  and  Then,  there  exists  an  invertible  matrix  T 
such  that 

B^  =  TB^  ,  C^  =  C^T-^  (2.14a) 

0/(P/  -  T-^Vi^T)Rg^  -  0/(F/  -  r-iy/r)P/  =  O  ,  (2.14b) 

and 

(Ai  -  T-U^T)Rg^  =  {El  -  T-^E^T)Rg'^  =  0  (2.14c) 

0/(Ai  -  T-^A,T)  =  Og^iEi  -  T-^E,T)  =  0  ,  (2.14d) 

where  i?/  and  0^^  are  the  strong  reachability  and  observability  matrices  for  sys¬ 
tem  1. 


2.4.  Extendibility 

The  concept  of  extendibility  for  stationary  TPBVDSs  was  introduced  in  [10]. 

Definition  2.1:  The  stationary  TPBVDS  (2.1)-(2.3)  is  extendible  (or  input- 
output  extendible)  if  given  any  interval  [K,L]  containing  [0,V],  there  exists  a  sta¬ 
tionary  TPBVDS  over  this  larger  interval  with  the  same  dynamics  as  in  (2.1),  but 
with  new  boundary  matrices  Vi{K,L)  and  Vf{K,L)  such  that  the  weighting  pat¬ 
tern  1T^(A:)  of  the  original  system  is  the  restriction  of  the  weighting  pattern 
-kI^)  ^be  new  extended  system,  i.e., 

W^{k)  =  WK_]^{k)  for  1-N  <k  <N  . 


(2.15) 
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Theorem  2.3  [10]:  A  stationary  TPBVDS  is  extendible  if  and  only  if 


0,{Vi  -  ViE^E)R,  =0  (2.16a) 

0,[Vj  -  VjA^ A)R,  (2.16b) 


where  E^  and  A^  denote  the  Drazin  inverses  [13,  p.8]  of  E  and  A. 

From  conditions  (2.16),  by  using  the  E-,  A-,  E^-,  and  A^-  invariance  of 
Im(i2^)  [10]  and  the  generalized  Cayley-Hamilton  theorem  for  the  pencil  {E ,A} 
[9],  it  is  easy  to  check  that  for  an  extendible  stationary  TPBVDS,  the  weighting 
pattern  (2.12)  can  be  rewritten  as 


W{k)  =  1 


CViE^E^{AE^f-'^B 

CVfA^A^{EA^)-^B 


l<k  <N 
1-N  <k  <0  . 


(2.17) 


Given  an  extendible  stationary  TPBVDS  over  [0,V]  with  weighting  pattern 
VV)y(A;),  it  is  of  interest  to  ask  whether  it  is  possible  to  extend  this  TPBVDS  in  a 
consistent  way  over  intervals  of  increasing  lengths,  so  that  this  progressive  exten¬ 
sion  process  giv^  rise  to  a  unique  extended  weighting  pattern  VF(A:)  defined  for  all 
A; .  It  is  shown  in  [12]  that  for  an  interval  of  length  M  >  N,  if  we  select 


Vi  =  ViE^{E^)^  ,  Vf  =  VfA^{A^)^  ,  (2.18) 

as  the  new  boundary  matrices  over  this  larger  interval,  the  TPBVDS 
{C ,Vi,Vj  ,E ,A  ,B ,M)  is  an  extension  of  {C ,V^,Vf  ,E ,A  ,B ,N)  which  is  normal¬ 
ized,  stationary  and  extendible.  By  using  this  extension  procedure,  we  find  that 


W{k)  =  1 


C{ViE^)E^{AE^f-^B 

-C{I  -  {ViE^)]A^  {EA^)-’^  B 


k  >  Q 
k  <0  , 


(2.19) 


is  the  desired  extended  weighting  pattern. 


3.  Internal  Description  of  a  Weighting  Pattern 

The  matrix  V^E^  specifies  entirely  the  effect  of  the  boundary  conditions  on 
the  extended  weighting  pattern  ly’(A)  given  by  (2.19).  This  motivates  the  intro¬ 
duction  of  the  following  concept. 

Definition  3.1;  Let  {C ,E ,A  ,B ,N)  be  a  stationary  and  extendible 
TPBVDS.  Then,  P  is  a  projection  matrix  of  this  system  if 

0,PR.  =  0,{E«  Vi)R,  .  (3.1) 


The  extended  weighting  pattern  (2.19)  can  be  expresssed  in  terms  of  P  as 


W{k)  = 


I  CPE^{AE^f-'^B  k  >  0 
\-C{I-P)A^{EA^)-’^B  k  <0  . 


(3.2) 


Also,  by  using  (2.9),  (2.16),  (3.1),  and  the  fact  that  Im(P^)  and  Ker(Oj)  are  E- 
and  A  -  invariant,  it  is  easy  to  check  that  a  projection  matrix  P  satisfies 
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0,{PA  -  AP)R,  =  0,{PE  -  EP)R,  =  0  (3.3a) 

0,{P  -  PEE^)R,  =  0,({/  -  P)  -  {I-P)AA^\R,  =  0  .  (3.3b) 

As  is  clear  from  Definition  3.1,  one  particular  choice  of  projection  matrix  is 
P  =  F,  P^.  This  choice  is  not  unique  in  general.  If  P  is  a  projection  matrix,  so  is 
P  +  Q ,  where  Q  is  any  matrix  such  that  QRg  equals  zero. 

The  expression  (3.2)  for  the  extended  weighting  pattern  IV’(^)  motivates  the 
introduction  of  the  following  concept. 

Definition  3.2:  A  5-tuple  {C ,P ,E ,A  ,B)  is  said  to  be  an  internal  description 
of  the  acausal  weighting  pattern  W{k)  if  it  satisfies  (3.2)  and  (3.3),  and  if  {E ,A} 
is  in  standard  form.  Furthermore,  [C ,P ,E ,A  ,B)  is  minimal  is  it  has  the  smallest 
dimension  among  all  internal  descriptions  of  W{k). 

Given  an  acausal  weighting  pattern  W(/:),  a  possible  procedure  for  construct¬ 
ing  a  minimal,  extendible,  stationary  TPBVDS  {C ,V^,V^  ,E ,A  ,B ,N)  which 
admits  W{k)  as  extended  weighting  pattern  consists  therefore  in  dividing  the  real¬ 
ization  problem  into  two  steps.  First,  find  a  minimal  internal  description 
{C ,P ,E ,A  ,B)  of  W{k).  Next,  given  a  finite  interval  [0,iV],  find  some  appropriate 
boundary  matrices  F,-  and  such  that  the  corresponding  TPBVDS  is  extendible 
and  stationary,  and  such  that  P  is  a  projection  matrix  associated  to  these 
matrices.  The  following  result  guarantees  the  validity  of  this  two-step  realization 
approach. 

Theorem  3.1:  Let  [C ,P  ,E ,A  ,B)  be  an  internal  description  of  1F(A:).  Then, 
for  any  interval  length  N,  there  exists  matrices  F,-  and  F^  such  that  the 
TPBVDS  {C ,Vf,Vf  ,E ,A  ,B  ,N)  is  normalized,  extendible,  stationary,  and  has  for 
extended  weighting  pattern  1F(A:).  P  is  a  projection  matrix  of 
(C',F,-,  Fy  ,P,A  ,P,V).  Furthermore,  this  TPBVDS  is  minimal  if  and  only  if  the 
internal  description  {C,P,E , A, B)  of  W{k)  is  minimal. 

Proof:  Let 

F-  =  P{E^f  4-  crX{aE^  4-  A^)”^  (3.4a) 

Fy  =  (I  -  P)(A^)^  -h  X{aE^  4-  A^)-^  ,  (3.4b) 

where 

X  -  /  -  PEE^  -  {I  -  P)AA^  =  [I  -  P)EE^  4-  PAA^  -  EE^  AA^  ,  (3.5) 

and  where  a  is  any  scalar  such  that  aE^  4-  A-^  is  invertible.  The  relations  (3.4)- 
(3.5)  specify  a  TPBVDS  (C ,Vf,Vf  ,E ,A  ,B ,N).  By  direct  calculation,  it  is  easy  to 
check  that  F,-  and  Fy  are  normalized,  and  that  this  TPBVDS  is  stationary, 
extendible,  and  realizes  W{k).  Noting  that  O^XR^  =  0,  we  can  also  verify  that  P 
is  a  projection  matrix  of  the  TPBVDS.  By  construction,  the  minimality  of  the 
TPBVDS  {C ,V^,Vj  ,E ,A  ,B ,N)  is  equivalent  to  that  of  the  internal  description 
{C,P,E,A,B).  □ 

Given  an  internal  description  {C ,P ,E ,A  ,B)  of  the  weighting  pattern  W{k), 
the  following  result,  which  was  derived  in  [12],  shows  that  it  is  possible  to 
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characterize  the  minimality  of  this  internal  description  directly,  without  invoking 
minimality  conditions  for  an  associated  TPBVDS. 

Theorem  3.2:  The  internal  description  [C ,P ,E ,A  ,B)  of  W{k)  is  minimal  if 
and  only  if 


(a)  —  [Rg  P^s]  row  rank 


(b)  O. 


OsP 


has  full  column  rank 


(3.6a) 


(3.6b) 


(c)  KeviO,)  GMRs)  ■ 


(3.6c) 


By  analogy  with  the  weak  reachability  and  observability  matrices  which  were 
introduced  in  [9]- [10]  to  characterize  the  concepts  of  weak  reachability  and  obser¬ 
vability  for  a  TPBVDS  {C ,Vj,Vf  ,E ,A  ,B ,N),  the  matrices  R^  and  are  called 
the  weak  reachability  and  weak  observability  matrices  of  the  internal  description 
[C,P,E,A,B). 

It  is  also  shown  in  [12]  that  Theorem  3.2  implies  that  two  minimal  internal 
descriptions  of  a  weighting  pattern  can  be  related  as  follows. 

Corollary:  Consider  two  minimal  internal  descriptions  {C^  ,Pj,Ej,Aj,Bj), 
with  j  =  1,2,  of  the  same  weighting  pattern  W{k),  which  are  in  standard  form 
for  the  same  a.  and  /?.  Then,  there  exists  an  invertible  matrix  T  such  that  rela¬ 
tions  (2.14a),  (2.14c)-(2.14d),  and 

0}{P^-T-^P^T)R}=Q  (3.7) 

are  satisfied. 

The  procedure  that  we  develop  here  for  constructing  a  minimal  internal 
description  {C ,P ,E ,A  ,B)  of  W{k)  relies  on  the  introduction  of  a  new  transform, 
the  {s,t)  transform,  and  on  formulating  the  realization  problem  as  a  factorization 
problem  for  rational  matrices  in  s  and  t . 


4.  (s,i)-Transform  and  Rational  Matrix  Factorization 

One  difficulty  associated  with  the  use  of  the  2-transform  for  analyzing 
discrete-time  descriptor  systems,  is  that  since  the  dynamics  of  such  systems  are 
singular,  infinite  frequencies  cannot  be  handled  in  the  same  way  as  other  frequen¬ 
cies  [14].  This  motivates  the  introduction  of  the  transform 

CO 

H{s,t)  =  S  H{k)t^-ys’^  .  (4.1) 

k——oo 

It  can  be  expressed  in  terms  of  the  standard  2-transform  B[{z)  as 

H[s,t)=H{slt)lt  .  (4.2) 

From  this  observation,  we  see  that  when  H{s,t)  exists,  it  is  always  stricly  proper 
in  {s ,t)  in  the  sense  that 

lim  H{cs ,ct)  —  0  . 

C  — i-oo 


(4.3) 


8 


Note  however  that  it  is  not  necessarily  strictly  proper  in  s  and  i  separately,  so 
that  the  corresponding  z -transform  may  not  be  proper. 

In  the  following,  we  shall  restrict  our  attention  to  the  case  when  JI(z)  and 
H{s,t)  are  rational.  Then,  from  (4.2),  we  see  that  the  numerator  and  denominator 
polynomials  of  all  entries  of  H{s,t)  are  homogeneous.  Furthermore,  the  relative 
degree  in  s  and  t  of  all  entries  of  i.e.,  the  difference  between  the  denomi¬ 

nator  and  numerator  degrees,  is  exactly  one.  Thus,  the  transformation  (4.2)  has 
the  effect  of  transforming  rational  matrices  H{z),  proper  or  not,  into  strictly 
proper  homogeneous  rational  matrices  in  the  two  variables  s  and  t  with  relative 
degree  one. 


4.1.  Formulation  of  the  Realization  Problem 

In  the  causal  case,  the  z -transform  plays  an  important  role  in  the  solution  of 
the  minimal  realization  problem.  Specifically,  given  a  causal  weighting  pattern 
W{k),  the  minimal  realization  problem  is  equivalent  to  finding  matrices  {C,A,B) 
of  minimal  dimension  such  that  the  z-transform  W{z)  admits  the  the  factoriza¬ 
tion 

W{z)  =  C{zl  -  A)-^B  .  (4.4) 

For  the  case  of  acausal  weighting  patterns,  the  situation  is  more  complex.  If 
{C ,P ,E ,A  ,B)  is  an  internal  description  of  the  weighting  pattern  W{k),  and  if 
W^{k)  and  W^[k)  are  the  causal  and  anticausal  parts  of  W{k),  the  (s,i)- 
transforms  of  Wj  [k)  and  W^[k)  can  be  expressed  as 

OO 

Wf{s,t)=  ^CPE^iAE^f-'^Bt’^-^/s'^ 

k=i 

=  CPE^{sI  -  tAE^y^B  =  CP{sE  -  tA)-^B  (4.5a) 

Wb{s,t)=  f)  -C{I  -  P)A^{EA^fBy-ys’‘ 

h=—oo 

=  C{I  -  P)A^[sEA^  -  tiy^B  =  C [I  -  P)[sE  -  tAy^B  .  (4.5a) 


Note  that  Wf{s,t)  and  Wj(s,t)  do  not  have  in  general  the  same  regions  of  con¬ 
vergence.  However,  by  analytic  continuation,  it  is  possible  to  extend  their  domains 
of  definition  to  the  whole  plane,  while  using  the  same  notation.  This  yields  the 
factorizations 


Wf{s,t)  +  lFj(s,t)  =  C{sE  -  tAy^B  (4.6) 

[Wys,t)  Wb{s,t)]  =  C{sE  -tAy^lPB  {I  -P)B]  (4.7) 

{sE-tAy^B.  (4.8) 

Since  the  specification  of  an  acausal  weighting  pattern  W{k)  is  equivalent  to  the 


'wys,ty 

CP 

Wb{s,t) 

C{I  -  P) 
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specification  of  Wf  {s ,t)  and  Wj(s,i),  we  see  from  (4.6)-(4.8)  that  the  construction 
of  an  internal  description  {C ,P ,E ,A  ,B)  of  W{k)  involves  the  factorization  of 
three  homogeneous  rational  matrices  in  s  and  t,  instead  of  a  single  rational 
matrix  for  causal  systems. 

4.2.  Factorization  of  Homogeneous  Rational  Matrices  in  s  and  t 

The  above  discussion  motivates  the  following  minimal  factorization  problem: 
given  an  homogeneous  rational  matrix  function  H[s,t)  of  relative  degree  one,  find 
matrices  {K,D ,F ,G)  of  lowest  possible  dimension  such  that 

H{s,t)  =  K{sD  -  tF)-'^G  .  (4.9) 

An  important  feature  of  this  factorization  problem,  is  that  even  if  we  impose  the 
additional  requirement  that  {D,F}  should  be  in  standard  form,  i.e.,  that  there 
exists  a  and  /3  such  that 

aD  +  /3F  =  I  ,  (4.10) 

the  matrices  {K ,D ,F ,G)  are  not  unique.  To  insure  uniqueness,  a  and  jS  must  be 
chosen  a  priori.  They  can  be  chosen  arbitrarily,  as  long  as  is  defined. 

Theorem  4.1:  A  matrix  function  H{s,t)  admits  a  factorization  of  the  form 

(4.9)  if  and  only  if  it  is  homogeneous  in  s  and  t  with  relative  degree  one.  Under 

these  conditions,  if  {cx,/3)  is  a  pair  of  scalars  such  that  exists,  H{s,t) 

admits  a  unique  minimal  factorization,  up  to  a  similarity  transform,  satisfying 

(4.9) -(4.10).  The  dimension  r  of  this  minimal  realization,  i.e.,  the  size  of  D  and  F , 
is  given  by 

r  =  d{H{az,l  —  pz))  ,  (4.11) 

where  d(.)  denotes  the  usual  McMillan  degree,  and  where  H{az,l  —  Pz)  is  a 
strictly  proper  rational  matrix  in  z . 

Proof:  Necessity  is  obvious.  To  prove  sufficiency,  let  a  and  P  be  such  that 
H{a,—P)  exists.  Then,  consider  the  rational  matrix  H{az,l  —  Pz).  This  matrix  is 


strictly  proper  in  z 

because 

lim  H{az,l  —  Pz)  =  lim  H{a,—P)fz  =  0  . 

z — ►OO  2— »-oo 

(4.12) 

It  can  be  realized  as 

H{az,l  -  Pz)  =  K{zl  -  F)-^G  . 

(4.13) 

Now,  assume  that 

a  ^  0  (otherwise,  reverse  the  roles  of  D  and  F),  and  let 

w  =  a/ [at  +  ps)  ,  z  =  s /{at  +  Ps)  . 

(4.14) 

In  this  case 

s  —  az / w  ,  f  =  (1  —  Pz)/w  , 

(4.15) 

which  implies  that 

H{s,t)  =  wH{az,l  -  pz)  =  wK{zI  -  F)-^G  =  K{sD  -  tF)-^G  , 

(4.16) 
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where 

D=^{I-PF)/a.  (4.17) 

Since  there  is  a  one  to  one  correspondence  between  the  factorization  (4.13)  of 
H{az,l  —  jSz)  and  the  factorization  (4.16)-(4.17)  of  H{s,t),  the  dimension  and 
uniqueness  properties  of  these  two  factorizations  are  the  same.  This  implies  that 
minimal  factorizations  of  of  the  form  (4.9)-(4.10)  are  related  by  a  similar¬ 

ity  transform,  and  have  a  dimension  r  equal  to  the  McMillan  degree  of 
H{az,l  -  (3z).  □ 

Corollary:  The  factorization  (4.9)-(4.10)  is  minimal  if  and  only  if  [D  ,F  ,G) 
is  strongly  reachable  and  {K,D,F)  is  strongly  observable.  Furthermore,  the 
dimension  of  a  minimal  factorization  is  equal  to  the  rank  of  the  Hankel  matrix 
O^Rg,  where  and  are  the  strong  observability  and  reachability  matrices 
associated  respectively  to  {K,D,F)  and  [D  ,F ,G). 

One  unsatisfactory  aspect  of  Theorem  4.1  is  that  the  dimension  r  of  a 
minimal  factorization  of  is  characterized  in  terms  of  the  McMillan  degree 

of  the  1-D  rational  matrix  H{az,l  —  ^z),  and  not  directly  in  terms  of  It 

turns  out  that  it  is  possible  to  characterize  r  directly  from  H{s,t)  by  extending 
the  concept  of  McMillan  degree  as  follows. 

Definition  4.1:  Given  a  homogeneous  and  strictly  proper  rational  matrix 
H{s,t)  in  s  and  t,  the  McMillan  degree  of  H{s,t)  is  defined  as  the  degree  of  the 
least  common  multiple  of  the  denominators  of  all  minors  of 

It  was  shown  in  [12]  that: 

Theorem  4.2:  K  H{s,t)  is  factorizable,  i.e.,  if  it  is  homogeneous  of  relative 
degree  one,  the  dimension  of  a  minimal  factorization  of  H{s,t)  is  equal  to  its 
McMillan  degree. 


6.  Minimal  Realization 

In  Section  4,  it  was  shown  that  the  specification  of  an  internal  description 
{C ,P ,E ,A  ,B)  of  a  weighting  pattern  W{k)  yields  the  three  rational  matrix  fac¬ 
torizations  (4.6)-(4.8).  This  suggests  that  the  construction  of  a  minimal  internal 
description  of  W{k)  can  be  formulated  as  a  minimal  factorization  problem. 


5.1.  Dimension  of  a  Minimal  Realization 

Theorem  5.1:  The  dimension  n  of  a  minimal  internal  description  of  W[k) 
is  given  by 


n  =<H[Wf{s,t)  Wi{s,t)])  +  d{ 


W,(s,t) 


)-d{Wf{s,t)+  IVt(s,t)),(5.1) 


where  d(.)  is  the  generalized  McMillan  degree  introduced  in  Definition  4.1. 

Proof:  Let  (C ,P ,E ,A  ,B)  be  a  minimal  internal  description  of  W{k),  and  let 
OJ,  p  and  r  be  defined  as  follows 

oj=  d[C{sE  -tA)-^[PB  {I  -  P)B])  =  d{\Wf{sA)  kKj(s,0])  (5.2a) 
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d{ 


CP 

C{I-P) 


{sE  -  tA)~^B)  =  d{ 


Wf{s,t) 


) 


T  =  d{C{sE  -  tA)-^B)  =  d{Wf{s,t)  +  W,,{s,t)) 


(5.2b) 

{5.2c) 


F rom  tbe  corollary  of  Theorem  4.1,  it  follows  that  w,  p  and  r  are  the  ranks  of  the 
Hankel  matrices  O^R^,  and  O^R^,  respectively.  Then,  from  the  minimality 

conditions  (3.6a)-(3.6b),  R^  and  have  full  rank,  which  implies  that  to  and  p 
are  the  ranks  of  the  strong  observability  and  reachability  matrices  Og  and  Rg, 
respectively.  From  condition  (3.6c),  we  can  also  deduce  that  the  rank  of  OgRg 
equals  the  rank  of  Og  plus  that  of  Rg  minus  n ,  so  that 


r  —  a;  +  p  —  n  , 

which  implies  (5.1).  □ 

Example  5.1:  Consider  the  weighting  pattern 

k>l 


W{k) 


A:  <  1  ’ 


(5.3) 


(5.4) 


where  a  and  b  are  scalar  parameters  with  c  <  1.  From  Theorem  5.1,  we  find 
that  the  dimension  of  a  minimal  internal  description  of  ]F(A:)  is  given  by 


n 


ab 

s  —  at 


■])  +  d{ 


a 

s  —  at 
—ab 
s  —  at 


)  _  ^(il - ^) 

^  ^  s  -  at  ^ 


fl+l-l  =  l  for  b  =9^1 
|l  +  l-  0  =  2  for6=?il. 


When  b  ^  1,  3.  minimal  internal  description  of  W(k)  is 

C=al{l-b)  P=l/(l-6)  E  =1  A=c  5=1, 
and  for  6=1,  we  can  select 


1 

0 

fil 

c  =  [a  a]  P  = 

o 

_ 1 

o 

1 _ 

E  =  I  A  =al  5  = 

1 

rH 

_ j 

(5.6) 


(5.7) 


Thus,  6=1  can  be  viewed  as  a  singularity,  in  the  sense  that  the  dimension  of  a 
minimal  internal  description  of  W  is  2  only  when  6  is  exactly  equal  to  1.  □ 


5.2.  Minimal  Realization  Procedure 

One  interesting  aspect  of  Theorem  5.1  is  that  as  an  intermediate  step  in  the 
evaluation  of  the  dimension  n  of  a  minimal  internal  description  of  W{k),  we 
obtain  to  and  p,  which  are  respectively  the  ranks  of  the  strong  observability  and 
reachability  matrices  of  a  minimal  internal  description.  This  observation  was 
used  in  [12]  to  obtain  the  following  realization  procedure. 
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Step  1;  Construct  the  minimal  factorizations 

{Wf{s,t)  W,{s,t)]  C{sE  -  tA)-^[Bf  B^]  (5.8) 

Cf  ^  ^  ^ 

~  [sE  -tA)-^B  ,  (5.9) 

where  if  a  and  ^  are  such  that  Wj  (a,—/?)  and  VFj  (a,— /5)  are  defined,  the  pairs 
{E  ,A  }  and  {E  ,A  }  satisfy  the  normalization  condition  (2.4)  for  the  same  a  and  /?. 
Step  2:  Let 

B  =  Bj  A  Bjj  C  =  C j  A  Ci^  .  (5.10) 

From  (5.8)-(5.9),  we  find 

W[s,t)=Wf{s,t)AW^[s,t) 

=  C{sE  -  tA)-^B  =  C{sE  -  tA)-^B  ,  (5.11) 

so  that  {C,E,A,B)  and  {C,E,A,B)  are  two  factorizations,  in  general  non- 
minimal,  of  W{s,t).  The  minimality  of  factorizations  (5.8)  and  (5.9)  implies  that 
{C,E,A,B)  and  {C,E,A,B)  are  respectively  strongly  observable  and  strongly 
reachable.  By  decomposing  these  two  factorizations  in  strongly 
reachable/unreachable,  and  strongly  observable/unobservable  components,  respec¬ 
tively,  we  obtain 

_  Aj  Ag  _  Bi 

^  =  I  0  J,]  ^  [  0  J  (5-12) 

and 

El  E2  Ai  Ag  Bi 

C=10Cd  £=  [„  A=  5=  (5-13) 

Step  3:  From  (5.11),  we  find  that 

W{s,t)  =  C,{sE,  -  TA.r^B,  =  C^isE^  -  TAj-^Bg  ,  (5.14) 

where  the  factorizations  {C i,E i,A  i,B i)  and  {C 2,E ^,A ^,B 2)  are  both  strongly 
reachable  and  observable.  This  implies  that  they  must  be  related  by  a  similarity 
transformation,  i.e.,  there  exists  a  matrix  T  such  that 

Ci^CgT-^  E^=TE^T-^  A^  =  TA^T-'^  Bj  =  TBg  .  (5.15) 

The  matrix  T  is  given  by 

r  (5.16) 

where  and  denote  respectively  the  strong  reachability  matrices  of 

(Bl,Ai,Bi)  and  (B4,A4,Bg). 

Step  4:  The  matrices  C ,  E,  A  and  B  of  a  minimal  internal  description  can 
be  selected  as 


13 


E,T-^ 

* 

A,T-^  * 

E  = 

0 

E, 

E2 

A  = 

0  A  j  A2 

0 

0 

1 _ 

. . — "'1 

1 

0 

0 

_ 1 

C  =  [0  Cl  C2] 


B  = 


b 

0 


J 


(5.17a) 


(5.17b) 


where  *  indicates  an  arbitrary  block  entry.  The  role  of  the  similarity  transforma¬ 
tion  T  is  to  guarantee  that  the  component  which  is  common  to  factorizations 
(5.8)  and  (5.9)  is  expressed  in  the  same  coordinate  system.  Note  that  (5.17) 
corresponds  to  a  four  part  Kalman  decomposition  of  [C,E,AfB)  into  strongly 
reachable/unreachable  and  observable/unobservable  parts,  where  there  is  no 
unreachable  and  unobservable  component,  since  the  internal  description  that  we 
construct  must  be  minimal. 


Step  5:  The  matrix  P  is  obtained  by  solving 

O5  VR^  =  Hf  ,  (5.18a) 

where  Hf  denotes  the  Hankel  matrix  associated  to  the  causal  part  Wj{k)  of  the 
weighting  pattern,  and  setting 

p  =  .  (5.18b) 

Example  6.2:  Let 

|o  A;  =  1 

Then 


Wf  {s ,t)  — —t /s{s  —  t)  Wij{s ,t)  =  llis  —  t)  ,  (5.20) 

and  according  to  Theorem  5.1,  the  dimension  of  a  minimal  internal  description  of 
1T(A:)  is 


n=2+2-l=3 


(5.21) 


Since  oj  =  p  =  2,  we  can  also  conclude  that  the  minimal  internal  description  is 
neither  strongly  reachable  nor  strongly  observable.  To  obtain  a  minimal  descrip¬ 
tion,  the  first  step  is  to  perform  the  minimal  factorizations 


[Wf  Wj]  =  [. 


s[s  —  t)  s  —  t 


-  [1  l](sJ  -  t 


0  0 
0  1 


1 

-1 


-t 


Wf 

s(s  —  t) 

1 

1  1 

1  0 

-1 

IT, 

— 

— 

-1  0 

{si  —  t 

0  0 

1 

s  —  t 


0 

1 


(5.22a) 


(5.22b) 


which  satisfy  the  normalization  condition  (2.4)  with  a  =  1,  /3  =  0.  This  yields 
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B  = 


1 

0 


a  =  [1 0] . 


In  this  case,  we  can  select  T  =1,  and 


G  =  [0  11]  E  =  I  A  = 

- 1 

O 

o  o 

o  % 

_ 1 

B  = 

-1 

1 

0  0  ij 

0 

(5.23) 


(5.24) 


where  *  denotes  an  arbitrary  entry.  P  is  obtained  from  (5.18),  which  yields 


*  *  * 


P  = 


0  1* 
10* 


(5.25) 
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